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Abstract: In the context of brane gas cosmology in superstring theory, we show why 
it is impossible to simultaneously stabilize the dilaton and the radion with a general gas 
of strings (including massless modes) and D-branes. Although this requires invoking a 
different mechanism to stabilize these moduli fields, we find that the brane gas can still 
play a crucial role in the early universe in assisting moduli stabilization. We show that a 
modest energy density of specific types of brane gas can solve the overshoot problem that 
typically afflicts potentials arising from gaugino condensation. 
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1. Introduction 

A major success of brane gas cosmology (BGC) is the utilization of stringy effects to explain 
the origin of the hierarchy of dimensions. In the seminal proposal of Brandenberger and 
Vafa ly], it was argued that in the early universe all directions could fluctuate about the 
self-dual radius due to the presence of both winding and momentum modes. The argument 
asserts that strings will generically intersect in (3+l)-dimensional subspaces, so that such a 
subspace will lose its winding and subsequently expand into the large directions we observe 
today. This scenario was mathematically realized by Tseytlin and Vafa in the context of 
dilaton gravity ||, and has since been extended to include the effects of a gas of -Dp-branes, 
where Alexander, Easson, and Brandenberger || argued that such a gas would result in a 
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hierarchy of extra dimensions. Namely the original 9-dimensional spatial manifold should 
decompactify into a hierarchical product space of x T2 x T%. 

Subsequent investigations indicate that wound strings provide a mechanism for isotropiza- 
tion Q and stabilization || of the the compact dimensions, and that the mechanism works 
on toroidal orbifolds || . The framework for these results is usually the low energy effective 
action of type IIA string theory, where the salient differences from general relativity are 
a massless dilaton and the dynamics of extra dimensions. These differences lead to the 
result that negative pressure in the compact directions, due to wound strings, results in 
contraction, not acceleration. The dilaton is assumed to have no potential other than that 
which is induced by its coupling to the bulk string frame Lagrangian e~ 2< ^{R + (V</>) 2 ) and 
possible D-brane sources; most successes of BGC rely on the dynamical running of the 
dilaton toward weak string coupling, j s <1. 

On the other hand one would like to stabilize the dilaton at a value where g s is still 
large enough to be consistent with gauge coupling unification |7[]. Moreover if g s becomes 
too small, the interactions between strings become too weak to allow the annihilation of 
winding modes in three dimensions, where the space should be allowed to grow ||. (See 
also || for a discussion similar issues.) Rather there is only a window of finely-tuned initial 
conditions consistent with three dimensions ultimately growing to be large. A third reason 
that the dilaton must not continue to roll to arbitrarily small values is the constraint from 
fifth force experiments and null searches for time variation of physical constants [1C]; these 
preclude the dilaton from continuing to evolve at late times. 

For these reasons it is imperative to reconcile brane gas cosmology with the stabilization 
of both the radion and the dilaton. In [11] it was shown that using just the string winding 
and momentum modes this is not possible. We therefore first investigate whether by 
including more general string and brane states one can achieve such a stabilization. (A 
similar but less general analysis has been done in fl2| .) 

Our result, in the context of superstring theories, is that a general gas of D-branes and 
strings cannot stabilize both moduli, although they can stabilize one linear combination of 
them. However, the brane gas can still play an important role in the process of stabilizing 
both fields, due to the overshoot problem fl3|1 . A much-studied mechanism for stabilizing 
the radion involves adding racetrack potentials coming from gaugino condensation (and 
possibly an antibrane [0]). The Minkowski minimum of these potentials is typically sepa- 
rated from a runaway (decompactification) direction by a very small barrier, which would 
always be overcome by the inertia of the fields if their initial conditions were not finely 
tuned to be close to the desired minimum. One of our main observations is that a gas of 
brane winding modes can very robustly solve this problem by slowing down the modulus 
as it rolls down its steep potential. 

Our plan is as follows: in Section 2.1 we describe the BGC scenario and motivate 
the dimensional reduction procedure to obtain a d-dimensional theory of gravity with two 
scalar fields (the dilaton and radion), with an effective potential coming from the brane 
gas. In Section ^3] we discuss some of the features of the effective potential; namely, we 
show that provided the dilaton is stabilized by some other mechanism, branes can stabilize 
all the extra dimensions. Section |3| presents our no-go theorem showing that under the 
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given assumptions, there exists an unstabilized direction in the moduli space of the dilaton 
and radion no matter what modes are included in the gas of D-branes and strings. In 
particular we also show that the presence of massless F-string modes do not help in lifting 
the runaway direction. In section || we consider the combined effect of the brane gas with 
a superpotential, such as would arise from gaugino condensation and antibranes, and show 
that the brane gas can provide a remedy for the overshoot problem. We give our conclusions 
in section |E|. Technical details are given in the appendices. 



2. Effective Brane Gas Cosmology 

2.1 Supergravity coupled to Strings and Brane Sources 

A starting point for BGC is type IIA string theory compactified on a 9-dimensional toroidal 
background, which may be thought of as the result of compactifying M-theory on S 1 . The 
low-energy bulk effective action of this theory is given by 



1 



S "° = 2^j d " x ^ Z O e ' 2 * [ R + 40 MA, V,#V K * - j^H fm H^"\ , (2.1) 

where G is the determinant of the ten-dimensional background metric G^ u , <p is the dilaton, 
H is the field strength corresponding to the bulk antisymmetric tensor field B^, and k is 
the D-dimensional Newton's constant. For simplicity we ignore any flux contributions, and 
take H = 0. We envision this analysis to apply in the late-time era of BGC [l], ^, ||, 15], 



an epoch where the extra, compact dimensions are expected to be isotropized Ej, and 
winding modes in the large directions have annihilated. Thus, we consider a spacetime 
consisting of a flat, d-dimensional FRW universe, and an isotropic compact subspace of n 
extra dimensions 

ds 2 = G MN dX M dX N = g^dxV + b 2 (t) lmn dy m dy n (2.2) 
= -dt 2 + a 2 (t) dxidx 1 + b 2 (t) dy m dy m , i E {1, . . . , d}, m e {1, . . . , n} (2.3) 

where y m are the coordinates of the n extra dimensions. The total action comprises the 
above bulk action (2.1) and the action of all matter present. Sources are included by 



adding matter terms for both the strings (p s ) and Dp-branes (p p ). Owing to the different 
world-sheet couplings between the dilaton and the branes and strings, the matter action 
has the form 



S m = - d D xV^G p s + (2.4) 




T MN = — ^j^j • (2.5) 

We continue the construction of late-time BGC by considering separate species of 
strings and branes, each possibly having excited momentum (in the case of branes also 
known as "vibrational modes" [ 16 1 ) in the large or compact subspaces, but having winding 
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modes only along the compact directions. Then one can show (see appendix ^) that the 
stress energy tensor for the strings and branes simplifies to 



Ps + e 'Ppy 
P8 a b 



b 



i 
i 

J2&i ^ Pi e- ai<l> a~ d{1+iJi) b- n(1+Cli) 



0), 



(2.6) 
(2.7) 



In the preceding expressions the summation is performed over the relevant modes con- 
tributing to the gas of strings and branes, P and p being the sum-total pressure along 
the large and compact directions respectively. «j = for string sources, an = 1 for brane 
sources, and pi is the initial energy density for a particular mode, with effective equation 
of state pi = u>iPi, P{ = u>iPi. The values of u and to depend on the specific type of mode, 
dimensionality of the branes and the number of large and extra dimensions (see table [l]), 
but the important thing is that all the known modes can be described by these quantities. 

Variation of the action fl2.1| ) together with the matter action ( p.4| ) and metric ansatz (J2T 
results in the system of equations 




e^E 



\e*p 



(2.9) 
(2.10) 
(2.11) 

(2.12) 



l a(t) 



n 



(recall that d 



where we have introduced the shifted dilaton as cp = 2(f> — d- ( 

3 and n = 6), and a dot denotes differentiation with respect to time. Eqs. (pT£ - 2.12] ) 
are the string frame, or dilaton-gravity, equations of motion. Equation ( |2.9| ) is the 0-0 
Einstein equation; notice that in the string frame the kinetic term for the (shifted) dilaton 
contributes to the energy with apparently the wrong sign — this is due to the nonminimal 
coupling between the Ricci Scalar and the dilaton. The spatial components of the Einstein 
equations ( p.lC -2.11) show that the acceleration of the scale factor is proportional to the 
pressure, and thus the negative-pressure winding modes lead to contraction — this is the 
key ingredient of the Brandenberger-Vafa mechanism. Eq. ( |2.12 ) is the dilaton equation 
of motion. Equation ( |2.9[ ) is not dynamical, but is rather an equation of constraint, which 
can be used to determine the initial dilaton velocity to be 



v = ± 



e^E + d 



d(t) 
o(t) 



+ n 




(2.13) 
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It is customary to choose the negative solution since the string coupling g s = then evolves 
toward weak coupling, where a perturbative description is valid. Since all the terms under 
the square root are positive, the dilaton cannot bounce. 

The rolling of the dilaton, although important for the BV mechanism and the stabi- 
lization of the moduli fields, may also be deleterious to the BGC scenario. In ||, Easther, 
Greene, Jackson, and Kabat show that if the dilaton rolls too quickly, winding-mode an- 
nihilation may be suppressed, so that dynamical evolution leading to three large spatial 
dimensions is not favoured. The rolling of the dilaton also implies evolution of volume of 
the compact space. A conformal transformation on the metric may absorb the 4>-R cou- 
pling term, but this means the Einstein frame scale factors get additional time dependence 
from (f){t). This problem has typically been set aside (on the assumption that the dilaton 
will be stabilized at a later time) in discussions of stabilization of the extra dimensions || . 
However, in order to have a complete and consistent picture in the framework of brane 
gas cosmology one indeed needs to address the issue of stabilizing the dilaton along with 
radion stabilization, and this is what we devote the next few subsections to. 

The preceding observations also stress the utility of viewing gravity from the point of 
view of the four-dimensional Einstein frame, which is more intuitive than the 10D string 
frame. An effective four-dimensional action is achieved by conformally absorbing the dila- 
ton, integrating out the extra dimensions, and performing a second conformal transforma- 
tion to absorb the scale factor of the extra dimensions. The result is a minimally-coupled 
theory of BGC, where the original string and brane sources act as an effective potential 
for both the radion and dilaton fields. This approach was first advocated in ||i"7|| to study 
stabilization of extra dimensions in the presence of hydrodynamical fluids and was used 
to study string winding and momentum modes in jnj. We now generalize the analysis to 
include all possible string and brane sources. 

2.2 Effective Potential 

Upon performing dimensional reduction on both the string and brane sources, we obtain 
general relativity coupled to two scalar fields, the dilaton and radion, with an effective 
potential coming from the brane gas. As outlined in Appendix [A|, a string/brane source 
whose energy density behaves as p = p^-^Q-^+^^^-^+^i) ^ with equations of state uji 
and uji in the d large and n compact directions respectively, provides an effective potential 
in d + 1 dimensions 




fii = - ( -dui - nuji + 1 - a, ) \l d + ~ _ 1 ( 2 - 14 ) 

This is expressed in terms of the canonically normalized moduli ip and dilaton ip fields, and 
the Einstein-frame scale factor a of the d large directions, cc, parametrizes string (pn = 0) 
or brane (ck, = 1) contributions, and pi is the initial energy density of the ith component 
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of the brane gas. We work in Planck units, M pl 2 = 8ttGn = 1- The net effective potential 



will comprise several contributions of the form (2.14), depending on the type of excited 
modes; Appendix || discusses the equations of state, and the coefficients pi, vi for the 
various string and brane sources and the results are summarized in Table 1. 
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Table 1: A summary of the radion and dilaton couplings in the effective potential due to various 
species of gas. The spatial background consists of d large and n compact directions, with equations 
of states lo and Co respectively. 



2.3 Radion Stabilization 

To understand the effects of string and brane sources in late-time BGC, we now specialize 
to the case of three large directions (d = 3) with winding modes only in the compact 
dimensions. First suppose that brane sources are not present, so the effective potential for 
the system is given by contributions from strings alone — this emulates the setup of [jl], [5], 



11, [L5|, [TS|] . Three representative species of strings are considered, namely, W: strings with 
winding numbers in the compact direction (uj = 0, Co = — ^), Mq: momentum excitations 
in the compact directions (lo = 0, lo = and M3: momentum in the large directions 
(00 = 4, Co = 0). Summing contributions ( |2.14| ), we obtain 

Vsia,^) = ^ e ( 1 -f)^ e V / 4^-3 + ^ 6e -(i + f)^-3 + pM ^-4 t (215) 

where p w , p Ms , p M parametrize the initial energy densities of the three kinds of compo- 
nents, B = , 2 .^ , and A = Let us assume that the dilaton has been stabilized by 

' n(n+2) ' n+2 •> 

an external potential and consider the effect of the string gases on the unstabilized radion. 
Taking the dilaton VEV to be (/> = and ignoring the M3 momentum modes, which anyway 
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gets redshifted by the expansion of the universe, the string gas effective potential ( 2.1o| ) 
becomes 



V.M) = a' 3 \ PN e&-V^ + pue-W^] . (2.16) 



Battefeld and Watson point out [11] that this is a stable potential for ip only if the number 
of extra dimensions is n = 1, in which case it reduces to V(a,ip = 0) ~ Jy. This can be 
considered a source of dark matter, similar to the string-inspired example of Gubser and 
Peebles fll9| . However, in the case of n > 2, jn]] points out that the effective potential 
behaves as V(a,ifi) ~ e~ a ^/a 3 , so that the radion also runs away to oo. Since n = 6, one 
sees that the presence of strings cannot stabilize the dilaton or the radion. 

We note that in [20, ^] massless string states were invoked to obtain stabilization of 



the moduli. However, the former are not present in the type II string (being removed by 
the GSO projection). Although, they are present at the self-dual radius in the heterotic 
string, 1 additionally, [El] requires quantized modes of the D-string to achieve complete 
stabilization of all moduli. It is not clear to us that the D-string can be quantized in the 
same way as the fundamental string. 

Let us therefore consider whether extending the analysis of |y|] to the case of general 
brane sources can solve the problem of moduli stabilization. Consider the contributions 
to V e s coming from p-branes wrapping the compact dimensions (u = 0, u = — —, denoted 
N), and momentum modes in the compact dimensions (u> = 0, u) = — , denoted M). The 
net effective potential from equation ( |2.14j ) is 



V p (a,<p,ip) 



)v / B^>-f )>]& ■ a - e -(l+f )v / BV p "( 1 +f ) V / 4^ 



(2.17) 



This scenario is similar to those analyzed in g, ^, [l6|, p!q] . As we will now explore, 
the result of including a gas of branes is the improved stability of the radion. Inspection 
of ( |2.17| ) reveals that provided p > Q, all internal directions will be stabilized, since there 
are both rising and falling exponentials depending on tp: 

V p (a, = ^ e <P-3)V5*o- 3 + p^e-^^a- 3 . (2.18) 

where again we have assumed = 0. This has a nontrivial minimum close to ip = 
provided that p > ^. Since string theory requires n = 6, the presence of (p > 3)-branes in 
the compact directions will stabilize the moduli. 

However a more detailed analysis may be necessary to realize these stability conditions: 
According to the heuristic argument of Alexander, Brandenberger, and Easson ||, winding 
modes will generically intersect in 2p+ 1 dimensions, so that only objects with p < 2 should 
remain wound in the 6 compact directions. In this case, the stability requirement will not 
be satisfied. On the other hand, a quantitative investigation should account for the larger 
phase space once the a directions have grown large, thus decreasing the probability of 
annihilation, and perhaps leaving some extended objects with p > ^. As well, such an 



x We thank Subodh Patil for discussions on this point. 
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analysis should be carried out within the product space x T2 x T% argued for by Q, not 
the T n x 73 topology we have considered. We do note though that, as long as the shape 
moduli is frozen, the exact shape of the tori does not matter and our stability analysis for 
the volume still applies. 



3. No-Go Result in Type II String Theory 

In the previous section we saw that the radion can be stabilized by a brane gas when the 
dilaton is assumed to be fixed; similarly one can show the corresponding result when the 
roles of the dilaton and radion are interchanged. One may naturally wonder whether both 
of these moduli can be simultaneously stabilized using the most general combination of 
string and brane sources. As we now show, this is impossible to do with the conventional 
(winding, momentum or oscillator) string and brane excitations. The argument is made in 
two steps, starting first with gases where each string or brane has only one kind of excitation 
("simple states"), although different species of strings or branes are allowed to co-exist. We 
then extend the argument to the more general case where individual components of the 
gas have more than one kind of excitation ("mixed states"). 

3.1 "Simple States" 

We consider the situation when the strings/branes have nontrivial wrapping of only some 
extra dimensions, i.e. it doesn't wrap the large dimensions. They thus appear point-like 
to 4D observers and redshift like nonrelativistic dust, a -3 , corresponding to d = 3, u = 
in ( [2.14 ), The equations of motion for the radion and the dilaton in the presence of such 



sources are 

f + 3H<p = -™*&& (3.1) 

^ + 3i^ = -^#^ (3.2) 
dip 

with 

V eS (^ VO = a- 3 Pi e 2 ^+ 2 ^ (3.3) 

i 

where the sum runs over all possible string and brane states. As summarized in Appendix |B], 
the exponents m, Vi depend only on the effective equation of state parameter u along the 
extra dimensions ( p. 14 ), and the coupling exponent a of these states to the dilaton in the 



string frame: 

5gas,i = - I d D x ^e- a '* Pi b- n(1+ ^ a~ 3 . (3.4) 



where a = 0, for the fundamental strings since the Nambu-Goto action does not contain 
any dilaton coupling, while for branes a = 1, originating from the dilaton coupling in the 
DBI action. 

After performing the conformal transformations involving the radion and the dilaton 
(see Appendix [A|) the above action gives rise to the effective potential (2.14) in the 4D 
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Einstein frame with 

1 

Mi 



2^+2 V y 2 J ' V 2(n + 2) V 2 

As noted earlier, the value of u> depends upon whether the mode in question has 
winding, momentum or string oscillations. To analyze the stability of a potential which 
is a sum over such modes, we will use the technique of |2^]: we identify the directions in 
the ip-ip plane in which there exists a rising exponential contribution. If such directions 
are sufficiently numerous, the system is completely stabilized. An exponential of the form 
p^mip+2viii> r j seg mos t steeply along the direction cos9i(p + sm9iifj where tan#j = Vi/ Hi. 
In the range 

9 = {0i- tt/2, 6i + vr/2) (3.6) 

there is a rising potential (wall) while along the other half-plane the potential asymptoti- 
cally falls to zero. Since our potential is a sum of exponentials, it is clear that: 

(I) There can be at most a single local minimum and no local maxima. 

(II) Such a minimum exists only if the potential grows in any direction away from the 
minimum. Thus there must exist angles 9i for which the ranges of angles in ( |3.6| ) cover the 
entire plane. 

By looking at the different directions of steepest ascent of the exponentials it is easy to 
verify whether (II) is satisfied. Curiously, for brane sources (cq = 1) we find a result that 
is specific to d = 3 large dimensions: the direction of steepest ascent is the same (modulo 
7r) for winding, momentum or any other modes. For the winding modes, it is given by 

t,n9 = ^ = f- =► 9=* + i°> P>3 (3.7) 
Hi V n 6 I 7T, p < 3 

while for momentum modes 

9 = ~ (3.8) 

for all p, as illustrated in figure 0. Thus using, say, a gas of 2- and 6-branes in type IIA 
theory, or a gas of 1 and 5-branes in IIB theory, one could stabilize all directions of the ip-ip 
plane except for those orthogonal to the direction of steepest ascents, that is, 9 = | ± | 
(see figure). Along these directions the potential is flat, and there is a zero mode. 

To lift the flat direction, one has to incorporate fundamental string sources, with 
(Xi = 0. According to Table p] for F-strings, the angle 

taxi = H = Jl^i±i (3.9) 

Hi V n Cji + i 

depends on d)j, and so the analysis has to be done separately for each case. For the string 
momentum, oscillatory and winding modes, & = 1/n, and — 1/n respectively. From ( |3.5| ) 
one then finds 

$mom = 77 ) ^osc = 77) ^wind = 77 (3.10) 

2 6 o 

All of these modes provide an ascending potential in the direction of 9 = — ^, which is the 
flat direction when only brane sources are present, but not along the opposite direction 
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runaway 
NS5-brane 




F-wind 



F-mom 



F-osc 



Figure 1: Directions of steepest ascent in the cf>-tfj plane for contributions from different brane gas 
sources, described in the text. Long dashed lines are for the hypothetical F-string sources with 



u) = ±1. All marked angles are 30°. Short dashed lines are for the sources described in section 3.3 



(see figure |l]). Thus after including the string modes, there is no longer a zero mode, but 
there is a runaway direction along 6 run = 

In fact we can make an even stronger statement. Suppose that there are other string 
sources we may be unaware of; nevertheless their equations of state should satisfy the weak 
energy condition — 1 < & < 1. Using the Table 1 entries for general string sources and 
varying Cj over this range gives an angle of steepest ascent in the range 

tt _ tan" 1 — < 9 < tan" 1 — -124° < 9 < 19° (3.11) 

5 - - 7 ~ v ) 

which again fails to lift the 6 run direction, as shown in figure |l]. 

An alternative understanding of the moduli instability can be directly inferred by a 
reparametrization of the effective potential (2.14) in terms of new fields \ = \[B^) + \ j ^ip 



and 7] = y -jip — \/Bip. The result for an arbitrary string and brane gas with three large 
directions and uj = is 

V eS ,i(a = 0) = Pi eH^+D+D* e -(f a' 3 (string) (3.12) 
V eSt i(a = 1) = Pi e -"(^+Dx a" 3 (brane) (3.13) 

Through a combination of sources it is possible to to stabilize the x mode; however string 
sources will only cause ij to grow, and a brane gas does not couple to r/. Thus 17 = 
4 if) — ^f~B(p is the unstable direction in field space, in terms of the unit vectors tp, ip. 
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This direction corresponds to the line 




(3.14) 



which coincides precisely with the principal runaway direction identified previously. 

The above result is consistent with ref. [^] , which used a perturbation analysis in the 
BGC scenario to show that the inclusion of branes alone is not enough to stabilize both 
the dilaton and moduli fields. Thus some other potential is needed to stabilize one of the 
moduli. Given such a potential, BGC does provide a mechanism of stabilizing the other 
degree of freedom provided that branes with p > § are present. 

The factorization of the effective potential is a coincidence of having d = 3 large 
dimensions, as can be seen from the nontrivial dependence on d in eq. ( |2.14j ). For scenarios 
other than d = 3, the gas of strings and branes is able to stabilize both fields. 

Finally, we note that in the above analysis we did not consider branes or strings which 
wrap some of the large three dimensions (u < 0) ; these do not give any additional leverage 
for stabilizing the radion. 

3.2 "Mixed States" and Massless Modes 

So far we have only considered states which are purely oscillatory, winding or momentum 
modes. More generally, strings could have a combination of such excitations. Can such 
mixed modes help in stabilizing the moduli? The answer, unfortunately is no. For massive 
modes the reasons are similar to the case of the simple states. The massless modes 2 have 
to be analyzed separately but they do not alter the conclusions. 
In the string frame, the string spectrum is 

2 

m F = 9" + NoscM * + ( 3 - 15 ) 

where M s is the string scale and the three terms on the right hand side correspond to the 
momentum, oscillatory and winding pieces respectively. The source action for the strings 
(|3.4) then becomes 

Sstr = / d D x ^a- 3 b~ n ^/m 2 F (b)+p 2 (3.16) 

where p is the momentum along the non-compact directions. After performing the dimen- 
sional reduction and conformal redefinitions, as usual we find that the effective potential 
for the canonical radion and dilaton coming from a gas of such states is given by 

V ef f(ip,^)=p = nE(ip,^) (3.17) 

where n is the number density, and E(tp, <j)) is the energy of these states which depends on 
both the dilaton and the radion. Since we already know how the exponents look like for 

2 Although, in the type II theory these are excitations of the tachyon which are removed by the GSO 
projection, such massless winding states are allowed in the heterotic string theory. 
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individual momentum, winding and oscillatory modes, it is easy to see that for the more 
general case the energy is just given by 



E(ip,(p) = ^Aff [m 2 e 4 0* TO ¥>+«W) + w 2 e 4(nw<p+v v ,il>) + N osc e*(v°v+ v °^)] +p 2 = ^m 2 F +p 2 

(3.18) 

with 

1 , 1 

Mm = ; Ho = and fi w = 

and 

u m = ; u = ^ and ^ = (3.19) 

First let us focus on the massive modes, for which one can ignore the momentum 
p. The key observation is that since N osc > for massive modes, the effective potential 
obtained in the Einstein frame must still satisfy conditions (I) and (II) above in order to 
have a local minimum. Since again such potentials can have at most one minimum and no 
maximum, if there exists a minimum, the potential has to keep rising along any direction 
as one tends towards infinity. Thus to determine whether there is a minimum, it suffices to 
investigate the behaviour at infinity in the cp-i/j plane. Going far enough toward oo along 



a generic direction, one of the three terms in ( 3.15 ) will dominate, and then our previous 
analysis applies, which assumed the presence of only one term in a given source. In the 
special direction where b remains constant, no one term dominates, but they all remain 
proportional to each other, behaving like a single term, so again the previous analysis 
remains valid. 

Next let us focus on the massless modes, for example, the ones considered in pj] . In 
this case N osc < 0; and depending on the winding and momentum quantum numbers one 
could have ml ~ (m/b — wbM 2 ) 2 leading to a minimum at b 2 = (m/w)M~ 2 (which is the 
self-dual radius when m/w = 1) |21j| . For these modes one can easily verify that the mass 
function can be cast as 



m 2 F (ip,i(j) ~ e^' [e^ - e (3.20) 

where 

V/ = ^ + ^/> (3-21) 

(as one can find by carefully tracing back the conformal transformations) is really the string 
frame radion and 

<ff = \p- ^ (3-22) 

is the orthogonal direction. As one can see, the mass ( |3,20| ) and the potential have a 
minimum at tp' = and hence the massless states stabilize the ip' direction, as argued in 



1 21]. However, ip' also precisely coincides with the direction that could be fixed just with 
winding branes. Thus we are still left with the orthogonal runaway direction (ip' — > —oo) 
that we found earlier. 
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3.3 Exotic States 

We have seen so far that ordinary D-brane and string states are unable to stabilize both 
the radion and the dilaton simultaneously. We now briefly discuss how stabilization might 
be achieved using some less conventional kinds of branes. 



One kind of exotic state which has been considered [21| are quantized D-string modes. 
Whether it is justified to derive these from the Nambu-Goto action like for F-strings seems 
doubtful, since the D-string is a solitonic object, but for completeness we have derived the 
exponents corresponding to the different D-string modes 3 . Although oscillator excitations 
do provide a new direction in field space whose potential has a steep direction, this direction 
overlaps with ones from other more conventional sources, and do not affect our no-go result. 
The direction of steepest ascent for the D-string oscillator modes is derived in Appendix |C], 
and is shown in figure [l]. However if massless D-string modes are also allowed in the string 
theory spectrum, they can lift the runaway direction in conjunction with other modes, as 



has been argued in [21] 



Another possible source that provides an effective potential is the NS5-brane. Its 
tension behaves as T§ oc gj 2 , so an NS5-brane wrapping the internal manifold corresponds 
to a = 2, uj = =^-, uj = 0; this results in the coupling coefficients 



M f = _(£ - 2) - 1 and v£ = (10 - n)J , 2 . (3.23) 
P5 y 2 J ^+2 5 V 7 yn(n + 2) v ' 

Thus, for the case n = 6, the potential rises maximally in the direction tp = —-^cf), 
corresponding to the angle 0F5 = 150° in figure pi This does not coincide with the runaway 
direction identified in figure ([!]), but is close enough so that NS5-branes in conjunction with 
strings or D-branes can stabilize all the moduli. 

4. Adding Superpotentials 

Since it is not possible to fully stabilize all moduli using the D-brane/string gas, we inves- 
tigate the dynamics of a system in which the brane gas is present simultaneously with an 
external stabilization mechanism. A typical potential which could arise in string-motivated 
supergravity theories is the one which is generated by gaugino condensation in an SU(N) 
gauge sector. Although it might seem redundant to consider partial modulus stabiliza- 
tion by a brane gas when there is already a potential at zero density, there could actually 
be several benefits: for example, the brane gas can prevent the problem of the moduli 
overshooting the desired minimum [13|, as we investigate in this section. 



4.1 Gaugino Condensation Potential 

We briefly review the derivation of the nonperturbative gaugino condensate potential in 
low-energy effective supergravity, starting with the 10 dimensional spacetime which is as- 
sumed to be a product of 4D noncompact external spacetime and a 6D compact internal 

3 See Appendix ^ for details. 

4 We thank Ali Kaya for pointing this out to us. 
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manifold. We limit our present discussion to the dynamics of the radion, ip(x). A similar 
discussion should apply for more than one moduli field, but for simplicity we assume that 
all other moduli (i.e., complex structure and dilaton) have been stabilized. The radion 
appears in the full metric as 

ds 2 = g^dx»dx v + e 2 ^g mn dy m dy n (4.1) 

In supergravity, the radion is the real part of a chiral field, 

T = X + %Y = e 4 ^ + iY (4.2) 

Dimensional reduction of the supergravity action yields an effective four dimensional theory 
of gravity coupled to the complex scalar field T(x) 



S = M 2 I (/'.r, 



j + K Tf d^T d^f - e K {lC TT D T W D T W - 3| W| 2 ) 



(4.3) 



where JC(T, T) and W(T, T) are the Kahler potential and superpotential respectively, while 
K T f is the Kahler metric given by 

JC TT = 7 ^L (4.4) 
TT dTdT y ' 

We have also performed a conformal transformation of the four dimensional metric: 

9w - e n V = ^g,u (4.5) 

The kinetic and potential terms for T are computed from K(T, T) and W(T, T), where 
the Kahler potential for T is K, = — 31n[T + T] while as in 113] we use the superpotential 

W = W + Ae~ aT (4.6) 

which would be obtained through gaugino condensation in a theory with a simple gauge 
group. For instance, for SU(N), a = 2tt/N. The constant term Wo represents the effective 
superpotential due to any fields that have been fixed already p3| , such as the dilaton and 
complex structure moduli. 5 

The scalar-tensor action then reads 

r p i 

(4.7) 

where the kinetic (K) and potential (V") terms are given by 





R 


1 d^xyf^g 


- + K-V 



8 T B^T 3 



5 Although recent authors have pointed out that a proper construction of the KKLT mechanism includes 
other non-perturbative contributions to the Kahler potential, we are primarily concerned with addressing 
the overshoot problem which still exists despite their inclusion {mI pfl. That is, the mechanism proposed 
still provides an attractor solution despite changes to the form of the Kahler potential. 
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and 

V = + 6X2 [ aA2 (aX + 3)e" 2aX + 3W Aae- aX cos(aY)] (4. 



To arrive at the potential fl4.8|) ) we have also included the potential energy coming from an 
anti-D3 brane (first term) as in [|l4|] , which is needed in order to have a nonnegative vacuum 
energy density. The coefficient E is a function of the tension of the brane T3 and of the 
warp factor, if there are warped throats |pq ] on the Calabi-Yau manifold. The exponent a 
is either a = 2 if the anti-D3 branes are sitting at the end of a warped throat. Otherwise 
a = 3 corresponding to the unwarped region. If a warped region exists, it is energetically 
preferred. 

The imaginary part of the Kahler modulus, the axion Y, has stable minima at Y = 
(2n + l)7r/a (assuming Wo-Aa > 0). We will integrate this field out and focus on the 
dynamics of the radion, whose kinetic term is 6 12M 2 (dtp) 2 , and whose potential becomes 

V = Ee- 4 "* - - {W Aae- aX - aA 2 e~ 2aX ) e'^ + \a 2 A 2 e' 2aX (4.9) 
2 (3 



It is convenient to rescale ip — > ipj\J2A so that the kinetic term is canonically normalized. 
The action becomes 

R (dip) 2 



S = M 2 I d?x^~g 



f Ee - 2a ^ - 1 (W Aae- aX - aA 2 e~ 2aX ) e~ 2a ^ + \ 2 A 2 e~ 2aX e~ 2a ^ 



(4.10) 



where 

«i = 4> a2 = ^ a3 = 7i' x = ^ ( 4 - n ) 

The potential has three distinct regions (see figure [2], solid curve). (1) For ip large and 
negative V(ip) is dominated by the antibrane contribution, if a = 3, or by a combination 
of the antibrane term and the term proportional to e" 2 " 2 ^ if a = 2. In either case, the 
potential is to a good approximation a pure exponential in ip, which will be relevant for the 
analytic solutions we discuss in the next subsection. (2) For ip ~ 1, the different terms in 
the potential are comparable, creating a minimum at V^min, followed by a potential barrier 
at VWx- (3) For ip ^> 1 the antibrane term again dominates, since the other terms are 
exponentially suppressed by e~ aX . 

4.2 Attractor solution with brane gases 

We now consider the effect of augmenting the vacuum potential in (4.10) with the con- 



tribution from a brane gas. Without the brane gas, the dynamics of the radion depend 
sensitively on the initial conditions. If we start with ip > f/Wx (the position of the bump 
in the potential), ip runs to infinity, where the extra dimensions are decompactified. 
Generically one might expect the radion to start closer to the Planck size with ip < 0, so 
that there is a possibility of reaching the stable minimum at ip = ip max . 

6 The kinetic part of the action for the radion can also be derived directly from the Einstein-Hilbert 
action S10 = / d 10 xy / —gR contained in the full 10D supergravity action. Using a consistent dimensional 
reduction ansatz for the metric of the form one obtains S4 = J d 4 x^—ge n ^[R — n(n — l)(dij)) 2 + ...]. 
A conformal transformation precisely of the form J4.5| ) is needed to convert the action S4 to Einstein frame, 
from which one recovers the kinetic piece of the radion given here. 
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g 

3 



> 



with brane gas 

vacuum potential 



Since the vacuum potential in the re- 
gion ij) < is well approximated by an ex- 
ponential, the radion quickly reaches the 
attractor solution discussed in p7J ; it tracks 
the minimum formed between the expo- 
nential potential and the rising part of the 
"brane-gas potential," shown as the dashed 
line in figure ||. This attractor behavior 
washes out the effect of initial conditions. 
As long as the attractor is reached be- 
fore the field has passed the position of 
the minimum, this will allow tp to settle 
into the minimum and avoid the overshoot 
problem. 

Let us recapitulate the details of the 
attractor solution. The rising part of the brane-gas potential originates from the winding 
modes of p-branes with p > 3. In this region the Friedmann equation and the equation of 
motion for if) read (in M p = 1 units) 

jl r*j 1 



Figure 2: Radion potential with vacuum gaugino 
condensate potential (solid line) and potential at 
nonzero brane gas density. 



±V> 2 + Ee~ 2a ^ + p p e 2u ^ 



with p p 



(4.12) 



(4.13) 



v p p p e 

respectively. The exponents v p for p-branes' coupling to the canonically normalized radion 
were derived in the previous section, 



2v n 



n(n + 2) 



n 

P ~2 




(4.14) 



where n = 6 is the number of extra dimensions. This kind of system was studied in [28], 
where it was shown that there exist tracking solutions in which the energy of the scalar 
field tracks that of the branes: 



E 

Pp 



«i(qi + v p ) - 3/4 N 





1 


~E 


i 








) 




— r 








.Pp _ 





(4.15) 



u p (u p + a\) + 3/8 _ 

This relation implies that both the potential and kinetic energy of the radion remain 
proportional to the energy density of the branes, 



r p p e p v 



(vp + ax, 



-K 



(4.16) 



3(1 +r) 

The steepest brane-induced effective potential occurs for the maximal value of p, p = 6; 
this provides the greatest resistance to expansion of the internal manifold and will be the 
most effective case for avoiding the overshoot problem. Eq. ( |4.16 ) also shows that for p = 6 
the ratio of kinetic to potential energy is minimized. For example, if a = 2 and p = 6 we 
have ai = y 7 ^, 2v 6 = y^S, leading to K/V = 12/23. 
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In passing we note that such tracking solutions correspond to a power law expansion 
of the universe 

t v (2/3)(l+^/ai) / t N 11/12 



«(*) = „„[ -I =„„(-] (4.17) 



The universe does not accelerate during this phase. However, as was found in |28|| , when 
the analysis is carried out including the dilaton, acceleration can be obtained. 

4.3 Addressing the Overshoot Problem 

The above discussion implies that the overshoot problem will be avoided in the presence of 
a brane gas so long as the attractor solution can be reached. This means that for a given 
initial value of ip, the initial energy density in the brane gas, p p e 2Up ^, must be sufficiently 
large. If not, the brane density is diluted too quickly by the expansion of the universe and 
the system evolves according to the vacuum potential. 

We have confirmed these expectations by numerically integrating the coupled system 
of Friedmann and radion equations, which we illustrate with a specific example. In the 
potential ( |4.8| ) we consider an antibrane in a warped throat, with a = 2. Its tension is 
tuned to give a Minkowski minimum as shown in figure [2], which illustrates the case where 
E = 0.00889, a = 2.1, A = 0.9, and W = 0.25. We first verified that indeed this potential 
suffers from an overshoot problem, shown in figure ||. Starting from an initial condition 
< — 0.17, the field runs away to do. 

Interestingly, overshooting can be prevented by initial brane densities which are many 
orders of magnitude smaller than the initial potential energy of the radion. Figure |] shows 
the evolution starting from exponentially large initial radion potential energy, with ipo = 
— 100 and p = 6, for several initial brane densities, parametrized by £ = p p e 2l,p ^ /Vo(ipo), 
where Vo is the potential of the radion alone, excluding the brane gas contribution. The 
result shows that even for initial brane gas energy densities which are only 10~ 18 ^0(^0)1 
overshoot can be prevented. For different initial values, the exponent log 10 (C) scales linearly 
with ipQ. This behavior can be understood analytically, as shown in Appendix D. The 
minimum required value of £ is given by 



log 10 C ~ -0.43 



(4.18) 



The intuitive explanation for this result is that the radion energy initially falls more quickly 
than that of the brane gas. What counts is not the initial ratio of brane gas to potential 
energy; rather it is the ratio at the time when if) is close to its nontrivial minimum. This 
mechanism has been pointed out in |2S!J (see also |3^, 31]) as a generic way of solving the 



overshoot problem, using general sources of energy density. Brane gas cosmology provides 
a concrete setting where this idea can be used advantageously. 

When the modulus has reached its stable minimum, we are still left with a gas of branes, 
whose energy density is comparable to the energy density in the scalar fields; otherwise 
the brane gas would not be effective in slowing the rolling of the modulus. At the bottom 
of its potential, the scalar field oscillates and and its energy density redshifts as a -3 just 
like the brane gas. The result is a matter dominated universe. We must assume that 
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20 40 60 80 100 120 

time 

Figure 3: Evolution of tp without brane gas, 
for several initial values tpo = —0.15, —0.16, 
• ■ ■ , —0.2, illustrating overshoot. 



20 40 60 80 100 120 

time 

Figure 4: Solutions with ipo — —100 and 
different initial densities of brane gas, near 
the borderline of overshooting. 



inflation begins some time after this in order to dilute the branes and reheat the universe. 
Work on smoothly connecting the modulus stabilization with the beginning of inflation is 
in progress. 



5. Conclusions 



In this paper we used dimensional reduction to derive the effective action for a gas of 
strings and p-branes, giving a contribution to the effective potential for the radion and 
dilaton. In a gas of strings only, this potential could stabilize the radion provided there 
was only one extra dimension, but not the dilaton. dilaton. Including p-branes allows 
for the stabilization of either the dilaton or radion if p > |. However, the brane gas is 
insufficient for stabilizing both moduli simultaneously, for the type II strings we consider, 
which have no massless winding modes. Rather, only a linear combination of the moduli 
can be stabilized by the brane gas. 

It thus seems likely that external potentials are needed for modulus stabilization. 
However the brane gas can still play an interesting role in helping the moduli settle into 
their typically shallow minima, avoiding the overshoot problem. An attractive feature of 
this mechanism is that the brane gas can initially be many orders of magnitude smaller in 
energy density than the potential energy of the moduli and still be effective in slowing the 
rolling of the moduli, since the brane gas energy redshifts more slowly. There is therefore 
no need for finely-tuned initial conditions. 

In this work we have ignored quantum corrections, as well as higher-derivative cor- 
rections to the dilaton gravity action. The first approximation is justified for weak string 
coupling, g s = « 1. In this regard, the runaway direction found in Section || corre- 
sponds to 4> —oo, showing that quantum corrections cannot lift this flat direction at 
large field values. Of course it is possible that such corrections could lead to a metastable 
minimum along the flat direction, which would be a loophole in our no-go result. 
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Note Added 



As this work was being finished, similar results were given in [32]. Our results were pre- 
sented the week before, at the McGill Workshop on String Gas Cosmology, 30 April 2005. 



A. Dimensional Reduction 



We provide here a brief review of the standard dimensional-reduction procedure, following 
the procedure of [JllJ, [£J. Our starting point is D-dimensional dilaton-gravity together 
with a generic contribution of gas. This system is described by 

f d D xV^Ge- 2 * (R + AG MN V M ^N^-^H^ ua H 



SlI = 



d D x, 



P 



(A.l) 
(A.2) 



for some initial density p^. The dimensional reduction procedure will focus on the string 
action ( A.l| ); but, by tracking the transformation rules, we can later also reduce the mat- 
ter components. We obtain an effective theory of BGC by first transforming the string 



action ( A.l ) to the Einstein frame through the conformal transformations 



Gmn — > Gmn = Q 2 Gmn, = e A ^ , A 

D -2A<p 



D-2 

A<t> I \: .\ r ,MX 
\ J ;MN 

(D-l)(D-A)(e~ A ^ 



R^R:R = e - 2A<p R - 2{D - l)e~ A(p [ e~ A(p ) G 



■M 



A<f> \ qMN 
N 



to obtain 



(A.3) 



S -> S 



<! L) .r \ ' U \ R - G MN V m4>V N 4>\ 



(A.4) 



where <j> is the canonically-normalized dilaton, and we have ignored flux contributions. 



We dimensionally reduce the action by integrating out the extra dimensions [11, 33]. To 
perform this last step we consider a string- frame metric of the form (|2.3|), split into d large 
directions described by and n compact directions described by 7 mn . For simplicity, 
we consider the geometry of the extra dimensions to be that of a torus, thus R[y m n\ = 
0. We use the following relations to isolate the scale-factor dependence on the extra- 



dimensions [11, 33. 34] 



R 



V-G = b n ^-~g 
R[G MN ] = R[9^} - 2n&- VV M V M 6 



n n 



(A.5) 
(A.6) 
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where, again, R[y mn ] = 0, n and b{x^) are the number and scale factor corresponding to the 
extra dimensions, and <jL„ is the metric of the non-compact directions. Since none of the 
terms in the action depend explicitly on the coordinates from the n extra dimensions, we 
integrate over these directions to get the low energy effective action of the d+l-dimensional 
theory 



Vn 



jd+l, 



-b n r u v^v^ 



(A.7) 



where V n = J d n y^j is the spatial volume of the n extra dimensions under unit scaling 
(6 = 1). 

A second conformal transformation and field redefinition of the action ( |A.7| ) is neces- 
sary to obtain the canonical form of the Einstein-Hilbert action. The conformal transfor- 
mation reuses the identities ( |A.3|) with 



9iiv — b g^ u — e v ' ^g^u > 



resulting in 



Seff = I d d+l x^[R[g^} - ^"V^V^ - ^"V^V^ 



(Ai 



(A.9) 



where B 



n (d+n-i) • Finally, the system is canonically normalized by identifying the 



4D Planck mass as = and by rescaling the fields as 



tP = M p ij , <p = M p 



(A.10) 



S, 



eff 



d d+l Xy 



Ml 



1 



1 



-fR[g^] - 2^v M ^v,v - -s^v^v^ ) . (A.ii) 



The net effect of these transformations is to rescale the scale factors and dilaton as 



-G 



a(t) -» a(t) = e d - lV M * e V ' ^a(t) 
b(t) -» b(t) = e^^k e~^~ 2 ^hb(t) 

<p{t) -> <p(t) = V2AM p( p(t), (A.12) 

Employing the above expressions, we may now express the contribution of a source behaving 
as 

p = /0 . a -d(i+" 1 ) & -™(i+^) ( A 13 ) 

in the D dimensional string frame, through the effective matter-action 



= I d d+1 xy/^gpie 



(A.14) 
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with M p = 1. The original theory of dilaton gravity together with string and brane sources 
can now be interpreted as a theory of Einstein gravity together with sources, plus two 
scalar fields corresponding to the dilaton (cp) and the moduli field (■0), this is the action 
of equation ( A.ll| ). As well, the source term (equation |A.2j ) now acts like an effective 
potential for the two scalar fields. The inclusion of different excited states will provide 
different effective potentials, and this freedom can be exploited in the search for a moduli- 
stabilizing potential. 



B. Equations of state 

In this section we derive the equations of state and the resultant coefficients for the brane- 
gas effective potential. Using the metric-ansatz ( |2.3[ ), we derive the gas pressure from the 
thermodynamic relation. 



P 

a ~ 5V 



(B.l) 

b=const. 



The volume is given by V = \/—G s = a d b n , while energy contributions are generically of 
the form E = a?b k = (a d ) 2 (&")", so that 

SV = b n 5(a d ) + a d 5(b n ) (B.2) 

SE = i^-5(a d ) + ^5(b n ) ( B - 3 ) 
a a a n b n 

j a?b k 5(a d ) E 

u- = up (B.4) 



P -_ 5 A 
a ~ 5V 

8E 

Pb -~6V 



-const 



a=const 



da d b n 5{a d ) V 
k aib k 5{b n ) ^E A r , 
■n^W) =UJ V =UJp ' (R5) 



IIUJ 



where we have made the identifications u = — 4 and uj = — — . Thus E = a^b k = a dw b 

a n 

The existence of winding and momentum modes for strings is a well-known result, and is 
the reason for the T-duality invariant spectrum of closed strings. Finding an embedding 
with quantized momentum modes of branes is less subtle because the T-dual of a wrapped 
brane results in a wrapped brane, not a momentum mode. However, we use the embedding 
described by Kaya fl6[| , which results in momentum modes in the compact direction with 
energy given by 

E « = wr (R6) 

where A n is an unknown eigenvalue for the n'th momentum mode (we choose A n > 0). The 
corresponding pressure due to this brane momentum-mode is 

P ° = W)' (R7) 

which is a positive quantity. 
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C. D-string oscillator modes 



For completeness we consider the naive quantization of D-strings, in case these modes could 
affect the no-go result for simultaneous stabilization of the dilaton and radion. Ignoring 
the 2-form gauge field that couples to the D-string, the spectrum of the D-strings looks 
identical to that of F-strings except for the replacement 

M s — ► M' s = e~^ 2 M s (C.l) 

The rescaling is again due to the dilaton coupling present in the DBI action for the D- 
strings. Provided we ignore N osc < modes, again it is sufficient to consider only the 
"pure" modes. A straight forward computation yields the following source actions 



Sd„ = J d D x y^a- 3 b-^ (C.2) 
S D ,osc = / d D x y^e^ 2 a- 3 b- n (C.3) 



and 

<S D ,wind = J d D x y/^fje-Pa-h-W . (C.4) 

The D-string momentum modes looks identical to those of the F-string momentum 
modes, and yield no new effect. The winding modes are the same as those obtained for 
Dl-branes, which have already been considered. The only qualitatively new contribution 
comes from the oscillatory D-string modes. Substituting a = 1/2 and Q = in ( |3.5| ) one 
finds 

7T 7T 

#D,osc = -77 - 77 (C5) 
2 D 

Again, this fails to stabilize the runaway direction. 
D. Solving the overshoot problem 

One can analytically estimate of what must be the initial ratio of energy densities in the 
brane gas and radion in order to solve the overshoot problem. If the initial energy density 
of branes is much smaller than the potential energy of the radion, the dynamical equations 
will be given by 

H 2 I ^ 2 + Ee~ 2a ^ (D.I) 

and 

+ 3Hip = 2aiEe' 2ai ^ (D.2) 
The radion rolls freely down the exponential potential and exact solutions are known |35j] : 



a 



t 1 ' 2 ^ and ~ ~ a 2ai (D.3) 



Thus the energy densities of the brane and the radion redshift in this non-tracking phase 
as 

Pp e 2v ^ ~ a -(3-^ P ai) whUe _ a -4a? (D.4) 
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Thus as long as 



3 — AvpCti < 4a 1 



(D.5) 



the brane energy density will catch up with the potential energy of the radion. We can 
calculate when this happens. The ratio of brane energy to radion potential energy is 



p p e 



2i> p ij) 



Vty) 



~ a 



4^ p ai+4af-3 



(D.6) 



and we want that this ratio to be 0(1), by the time the radion rolls to the minimum. Hence 
we need to start with an initial ratio such that 



c 



Ppoe 



where 



V(tp min ) 



ClQ 



,— 2aiipo 



« 



-4of 



(D.7) 



(D.f 



From (D.7) and <pM we find 



C = exp 



2a^o (l - ^p)] =► log loC « -0.43 [ 2 a^ (l - ^^)] 
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